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Abstract. The vacuum expectation value of the electromagnetic energy-momentum 
tensor between two parallel plates in spacetime dimensions D > 4 is calculated in the 
axial gauge. While the pressure between the plates agrees with the global Casimir 
force, the energy density is divergent at the plates and not compatible with the total 
energy which follows from the force. However, subtracting the divergent self-energies 
of the plates, the resulting energy is finite and consistent with the force. In analogy 
with the corresponding scalar case for spacetime dimensions D > 2, the divergent 
self-energy of a single plate can be related to the lack of conformal invariance of the 
\Q . electromagnetic Lagrangian for dimensions D > 4. 

o 

Two parallel, metallic plates separated by the distance L in vacuum, will interact 
o • due to the modifications of the quantum fluctuations of the electromagnetic field caused 
by the boundary conditions at the plates. The resulting force was first calculated by 
Casimir pQ who found it to be given by the attractive pressure P = — 7r 2 /240L 4 . Using 
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the conformal symmetry of the electromagnetic field in D = 4 spacetime dimensions, 
Brown and Maclay[2] later obtained the vacuum expectation values of all the components 
of the electromagnetic energy-momentum tensor 

where C = —(l/^F^ is the standard Lagrangian. While these expectation values were 
constant between the plates, the corresponding fluctuations of the separate electric and 
magnetic fields were found by Liitken and Ravndal to be in general non-constant and 
actually divergent as one approaches one of the plates [3]. These divergences are caused 
by imposing ideal boundary conditions valid for arbitrarily small wavelengths of the 
field. A physical boundary would only affect fluctuations down to a finite wavelength 
which is expected to result in an increasing, but finite value of the fluctuations near the 
plates. The quantitative effects of such more realistic boundary conditions have been 
investigated during the last few years but a complete and satisfactory description is still 
lacking [1]. 

Casimir forces in spacetimes with dimensions D > 4 were first systematically 
calculated by Ambj0rn and Wolfram [5J. For the electromagnetic field between two 
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parallel hyperplanes with separation L, the attractive pressure was found to be 

p.- (d - 1)(0 - 2 )QW) (2 ) 

where the factor D — 2 is the number of physical degrees of freedom in the field resulting 
from gauge invariance. If the energy density between the plates is constant, it would 
just be this pressure divided by the factor D — 1. This is the case when D = 4 and it 
is of interest to see if it holds also in the more general case D > 4. For this purpose 
we calculate in the following the separate fluctuations of the electric and magnetic 
components of the field which then allows us to find all the vacuum expectation values 
of the components of the energy- momentum tensor (CQ). 

Today these quantum effects could be of relevance for stacks of parallel branes 
where the electromagnetic field is replaced by one or more of the abelian Ramond- 
Ramond fields. Any divergent energy density would then have serious implications for 
the stability of such configurations due to the resulting large gravitational interactions. 

The electromagnetic field tensor = d^A u — d v A^ in D = d + 1 spacetime 
dimensions has d electric components Ei = F 0i and d(d — l)/2 magnetic components 
Bij = Fij. For the geometry under consideration, the simplest and most natural choice 
of gauge is the axial gauge n^A^ = where the unit .D-vector is normal to the 
plates. Taking this along the z-axis, we thus have A z = 0. The component A is no 
longer a free variable, but depends on the others via the Maxwell equation diF t0 = 0. 
It gives A = —A~ 1 diA i where the operator A = df. There are thus D — 2 independent 
degrees of freedom described by the spatial field components A{ where i ^ z. The full 
Lagrangian then follows as 

L = ^Jd d x [A< (<% - di A- 1 ^) Aj - A, (didj - 8ij A) A 3 ] (3) 

after a few partial integrations and neglecting surface terms. 

In order to quantize the system, we must solve the classical wave equation following 
from the Lagrangian. For this purpose we impose the boundary condition n^F^ v = 
at the plates. This is the same as for the MIT quark bag where it had a physical 
justification [6J. Here it is just taken for convenience. In the axial gauge it gives d z Ai = 
at the plates which is the Neumann boundary condition for each physical field component 
Ai(x) = Ai(t; x<r, z). We then have the general mode expansion 



which satisfies the wave equation and the boundary conditions. The factor y2/L is a 
normalization factor. In the mode sum we have dropped a term with n = since it will 
not contribute to any physical results after regularization. 

Quantization can now be done in the standard way. We introduce orthonormal 
polarization vectors eA normal to the wavevector and a longitudinal polarization 
vector along this direction. The coordinate components A in of the field are then 
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replaced by the polarization components (A\ n ,Az, n ). After quantization at t 
transverse components can then be written on the standard form as 



AxJkj 
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where oj\ = + k 2 z with k. 
the standard commutator 

[a An (k T ), ay n/ (k' T )] = Sx\'d nn >(27r) d ~ 1 S(k T - k^ 

However, the longitudinal component 



nn/L. The creation and annihilation operators now have 

(6) 
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contains an extra factor when the corresponding creation and annihilation operators 
have the same canonical commutator (JB]). The full field operator (J4]) is then expressed 
in terms of these new operators corresponding to definite polarization states. 

The field fluctuations between the two plates can now easily be calculated. As a 
simple example, consider E z = —d z A~ 1 djAj. If we isolate the mode with quantum 
numbers (n, k^), we find the operator 

[ikj)(-iuj n 
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acting on the vacuum state. Here we have used that A gives k^ + /c^ = lo\ in momentum 
space. We see that the transverse modes will not contribute here since they satisfy the 
orthogonality condition e A • k^ = 0. However, for the longitudinal mode we have instead 
bl • kr = kx and it will give a non-zero contribution. The derivative d z gives a factor 
k z and cos(mrz/L) — ► sin(n7r z / L) . For this mode alone we thus get the fluctuation 
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Including all the modes, we thus have for the full fluctuation of this electric field 
component 
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when we write k\ 
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For the other components we similarly find 

d d - l ki 



(27T 



u; n (<i 
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(10) 
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where there is an implied sum over the transverse index i. The magnetic field 
fluctuations can be obtained the same way and become 

d d ~ l k T \ kl , , „ J . o ^vt2; 
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(12) 
(13) 
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when we again sum over the indices i and j. We have also confirmed these results by 
performing the same calculations in Coulomb gauge instead of axial gauge. 

Using now a combination of dimensional and zeta-function regularization as 
previously used when D = 4 [7], we can write the result on the form 
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(14) 



{Att) d / 2 L d 

Here (rP) is the Riemann zeta-function while fjj(z/L) depends on the distance z from 
the plates. When the spacetime dimension D is even, it can be written on the compact 
form 

n / . \ D-l 



f D (z/L) 



71 



T(D) 



A. 



cot 9 (D = even) 



(15) 



where 9 = ttz/L. But when D is odd, no such closed expression is easily derived. 
However, using a different regularization based on the corresponding point-split Green's 
functions, one finds in general [S] 

1 



fn(z/L)= £ 



Ch(D,z/L) + Ch(D,1-z/L) 



(16) 



where Ch(D, z/L) is the Hurwitz zeta-function. When D is even, this can be shown to 
agree with (Tl5l) . 

The regularized fluctuations of the electric field normal to the plates thus become 



(El) 



(D - 2)V{D/2) 



(4tt) d / 2 L d 

while for the transverse components we find 

XD-2)Y{D/2) 



UP) - \f D {z/L) 



(Ef) = -2- 



(Att) d / 2 L d 

For the magnetic fluctuations we similarly have 

(D - 2)r(D/2) 



UP) + -f D (z/L) 



and 



Bf z ) = -(D-2] 



(B 2 ) = (D-3) 



(4tt)d/2 L d 
(D-2)T(D/2) 



Cr(D) - -f D {z/L) 



Cr(D) + -f D (z/L) 



(17) 



(18) 



(19) 



(20) 



i<3i \~ -J (An) D / 2 L D 

Notice again that in these expressions we have summed over the transverse indices i 
and j, each taking D — 2 different values. All these correlators are seen to diverge near 
the plates z —>■ or z — > L where the function f£>{z/L) diverges. This is the same 
phenomenon which has previously been seen in D = 4 dimensions |3J. 

The pressure betwen the plates due to these fluctuations is defined by P = (T zz ). 
From flU) we have T zz = B 2 7 — E 2 + C where now 



1 (D - 2)T(Dp) 

{ P ~ ~ 1) (47r) D/2 L D JD{Z/L) 



(21) 
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Together with the values for ( Ef) and ( Bf z ) from above, the ^-dependence from the 
function f D (z/L) cancels out in the pressure and gives the expected value ([2]). 

So far there are no inconsistencies in the obtained results. But when we now 
calculate the energy density £ = ( T 00 ) between the plates, with T 00 = Ef + E\ — C, we 
obtain 

8 = JD ^r/^ 2) [UD) ~ (D/2 - 2 )^/ L )] ( 22 ) 

The ^-dependence in the last term is non-zero when D > 4 and makes the energy density 
diverge like z~ D with distance z from the plates. As a result, the total energy of the 
system is infinite, a result which seems to be impossible to reconcile with the finite 
Casimir force ([2]). In fact, ([2]) corresponds to having a constant energy density equal 
to the first term in (1221) . This apparent inconsistency has been verified in a different 
approach based on Green's function methods j8]. 

It is tempting to explain this problem by the imposed boundary conditions. We 
have used the MIT boundary condition which is equivalent to letting the electromagnetic 
vector potential satisfy Neumann boundary conditions in the axial gauge. Had 
we instead imposed metallic boundary conditions, equivalent to Dirichlet boundary 
conditions for the vector potential in the axial gauge, the only change in the above 
results would be the replacement of the mode functions cos(mrz/L) with sm(nnz/L) in 
(jl]) so that — >• — /d in the above results. Needless to say, the problem would remain. 
Only for periodic boundary conditions, as for finite temperature, would the disturbing 
term be absent [9]. But this is not necessarily satisfying from a physical point of view. A 
more mathematical discussion of such divergences near confining boundaries has been 
initiated by Fulling but here only scalar fields are considered p]]. 

A physical explanation of the above conumdrum becomes apparent when we take 
the limit L — > oo and thus consider the quantum fluctuations around a single plate. 
From (1221) we then find the energy density 

e l = i D-2){D/2-2) 1 ^Q B (23) 

which is non-zero on both sides of the plate and diverges when we approach it. This 
situation is analogous to the diverging energy density surrounding a pointlike electron. 
It is intrinsic to a single plate and should not contribute to the interaction between the 
plates induced by the same vacuum fluctuations. To see the connection with the Casimir 
force, we should subtract the self-energy ff23l) for both plates from the full energy density 
( |22|) . taking into account both sides of each plate. We thus obtain the interaction energy 
density 

f (D/2-2)(L/(L-z)) D forz<0, 
£ = -( D - 2 ) (a\d/2td x \ ( R (D)-(D/2-2)f D (z/L) for < z < L, (24) 
1 ' { (D/2-2)(L/z) D for z>L, 

where now 
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f D (z/L) = Ch(D, 1 + z/L) + ( H (D, 2 - z/L) 



(25) 



It is seen to be finite everywhere, even at the plates. When integrating over the full 
volume, the ^-dependent terms cancels out as follows from 



Only the z-independent term in (1241) contributes and agrees perfectly with the total 
energy corresponding to the Casimir force. 

A similar and somewhat simpler system is the Casimir energy induced by a massless 
scalar field in the same geometry. One will then find a very similar result for the energy 
density as obtained here [8]. It diverges near the plates for all spacetime dimensions 
D > 2. Again this can be attributed to a divergent self-energy of each plate. However, 
when D = 2 there are no such divergences and zero self-energy. But this is also the 
dimension in which the scalar theory has conformal invariance. In higher dimensions 
D > 2 it is possible to make the scalar theory retain this invariance by adding a 
conformal term. The resulting, improved energy-momentum tensor [12] then contains 
an additional piece discovered by Huggins [T3] and makes it traceless. Including the 
Huggins term, the divergent part of the energy density corresponding to the last term 
in ( 1221) drops out as s first noticed by de Witt when D = 4[llJ. 

For the electromagnetic field we have used the canonical energy-momentum tensor 
dl]) which has the trace T M M = (4 — D)C. It is zero for D = 4 which reflects the 
well-known fact that Maxwell theory is then conformally invariant. There are then no 
diverences in the Casimir energy. Thus it is natural to relate the apparent inconsistency 
in the electromagnetic Casimir energy when D > 4 to the lack of conformal invariance. 
It does not seem to be possible to construct an improved energy-momentum tensor 
in this case because gauge invariance forbids the existence of any corresponding local 
Huggins term. From this point of view the divergent, electromagnetic self-energy can 
therefore not be removed. For this to be done, one needs a more realistic description of 
the boundary plates along the lines considered by others pE]. 

This work has been supported by the grants NFR 159637/ V30 and NFR 
151574/V30 from the Research Council of Norway. 
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